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ABSTRACT 

We present sensitive 850 fj,m SCUBA photometry of a statistically-complete sample of 
six of the most luminous (M B < -27.5; ^L„(1450A) > 1O 13 L ), z > 4 radio-quiet 
quasars, reaching noise levels (la ~ 1.5 mJy) comparable with the deep blank sky 
surveys. These observations probe the rest frame far infrared region (~150/xm), at 
luminosity levels for thermal sources comparable with those that IRAS permitted for 
low redshift quasars. One quasar (BR2237— 0607; z = 4.55) is detected at 850 /an with 
q ■ a flux of 5.0±l.lmJy (4.5<r), whilst a second (BR0019-1522; z=4.52) has a detection 

| at the 2cr level. When combined with our previous millimetre measurements of z>4 

Q\ , quasars, we find that there is a large range (5-10) in far infrared (FIR) luminosity 

(£fir) at fixed UV luminosity, and that the typical quasar has a Lfir and mass of cool 
(50 K) dust similar to that of the archetyepal low redshift (z=0.018) ultraluminous 
IRAS galaxy(ULIRG) Arp220 (L Fm - 5 x 1O 12 L ; M d (cool) - 1O 8 M ). If one 
assumes a fiducial FIR luminosity of 5 x 1O 12 L for for all quasars with Mb < —23, 
\ we find that ^15 per cent of the sources in the SCUBA deep surveys could be classical 

broad-lined radio-quiet AGN. Thus if one considers the observed ratio of Seyfert II to 
Seyfert I galaxies at low redshift and any contribution from totally optically obscured 
AGN, a significant fraction of the SCUBA source population will harbour AGN and 
hence the inferred star formation rates from submm fluxes may be overestimated if 
the active nuclei are bolometrically dominant or the IMF is top heavy. 
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1 INTRODUCTION 



Recent deep SCUBA surveys (5mail, Ivison fc Blain 1997 
Hughes et al. 1994 |Barger et al. 199S| ; [Ealcs et al. 199cj ) 
have revealed a significant population of discrete submil- 
limetre sources, amounting to ~1000 deg~ 2 above ~3 mJy at 
850 /jm. If these objects lie at redshifts greater than 1, they 
have a far infrared (FIR) luminosity Lfir ~ 10 12-13 Lq, as- 
suming a dust temperature of ~50K. Since the large negative 
.R'-correction of the dust thermal spectrum means that for 
1 z ^ 10, this luminosity is insensitive to source redshift 



(cf. Ar] |>220 in Figure 1). Prior to the SCUBA surveys, the 
only significant known high redshift population with com- 
parable bolometric luminosity were the quasars. However, 



the characterisation of the FIR properties of the high red- 
shift quasars is rather meagre. Following the first detections 
of rest frame FIR em ission from z>3 radio quiet quasars 
(McMahon et al. 1994), our current knowledge is still biased 



towards more lumi nous objects with 850 fim fluxes in the 



10-50 mJy range (Isaak et al. 1994 
Omont et al. 1996a|). 



McMahon et al. 1994 



The onset of nuclear activity seems to be closely con- 
nected with significant structure-formation events. At high 
redshift, black holes form in rare overden se regions and are 



fuelle d by their collapsing host objects (Haehnelt fc Ree, 1 



199S ) ; at low redshifts, intera ctions /mergers provoke in flow 



of gas to the galactic nuclei (Barnes & Hernquistl996), fu- 
elling central black holes. These processes would be expected 
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Table 1. Program objects, observing conditions and results. 



Object name 


z 


M B 


•Si. 4 GHz 


Rsa 


tint 


T 225 GHz 


T 850fj,m 


T45CVim 


-^850 


^450 








(3a) 




(mins) 


(mean) 






(mjy) 


(mjy) 


BR B0019-1522 


4.52 


-27.6 


<0.5 1 


<2 


108 


0.08 


0.21-0.26 


1.1-2.4 


3.3±1.6 


11±33 


PSS J0134+3307 


4.52 


-28.2 


<1.5 2 


<3 


108 


0.08 


0.27-0.34 


1.8-2.5 


0.4±1.5 


-3 ±54 


PSS J0248+1802 


4.43 


-27.9 


<1.5 2 


<5 


36 


0.10 


0.23-0.46 


2.4-3.2 


6.4±4.5 


-433 ± 264 


PSS J0747+4434 


4.42 


-28.2 


<0.5 3 


<5 


54 


0.08 


0.27-0.39 


1.8-2.5 


0.0±2.2 


72±105 


BR B1600+0729 


4.35 


-28.1 


<0.5 1 


<2 


108 


0.09 


0.28-0.33 


1.7-2.0 


1.8±1.6 


-7 ±55 


BR B2237-0607 


4.55 


-28.1 


<0.5 1 


<5 


144 


0.07 


0.29-0.41 


1.8-2.7 


5.0±1.1 


19±53 



Notes: Radio flux limits taken from (1) McMahon et al, in preparation; (2) Condon et al, 1998; (3) Becker et al ,1995 



to trigger a major starburst, though the relative prominence 
of either component, starburst or AGN, would depend on the 
precise dynamical details of the collapse or interaction. 

In this paper, we report SCUBA 850 fim photometry of 
a statistically complete sample of six of the most luminous, 
z > 4, radio-quiet quasars. Our previous work at 800 /im 
with the JCMT (Isaak et al. 1994) and 1.25 mm with the 
IRAM 30 m (McMahon et al. 1994, Omont et al. 1996a) and 
was only sensitive to relatively high luminosity sources. The 
aim of the current work is to carry out more sensitive obser- 
vations of a small sample of radio quiet quasars in an effort 
to make an initial estimate of the contribution of AGN at 
high redshift to the submillimetre source population. 

Unless stated otherwise, we assume Hq = hso x 
SOkms" 1 Mpc" 1 and q — 0.5 throughout this paper. 



2 OBSERVATIONS 

The quasars for study were selected from the APM sam- 
ples of high redshift quasars (Storrie-Lombardi et al. 1996, 
McMahon et al. in prep.) and PSS survey (Kennefick et 
al. 1995). From extant radio data as tabulated in Table 1, 
the quasars are radio quiet based on the definition, Log 
i?8.4<1.0, where Rsa=Lsa/Lb, is the ratio of radio and 
optical luminosities as defined by Hooper et al(1996), where 
spectral indices of —0.5 and —0.7 have been used to extrap- 
olate to the rest frame optical and radio fluxes. Moreover, 
based on the observed median spectral index of —0.7 found 
for radio loud quasars over the wavelength range, the ex- 
pected contamination at 850 is <0.1mJy is all cases. The 
most luminous z > 4 objects were selected in the HA range 
15-08 as determined by the telescope scheduling, providing 
effectively sparse sampling for Mb< — 27.5 

Observations were carried out during three 8 hr observ- 
ing shifts over the period 1998 August 16 to 19, with the Sub 



millimetre Comm on-User Bolometer Array (SCUBA) (Hol- 



land et al. 1998) on the 15 m James Clerk Maxwell Tele- 



scope. We used SCUBA in standard point-source photome- 
try mode at 450 and 850 /im. This involves placing the target 
on the central bolometer of each array, whilst 'jiggling' the 
secondary mirror in a 9-point pattern with 2 arcsec offsets, 
integrating for 1 second at each. Superimposed on this is a 
chop of the secondary mirror by 60 arcsec (in azimuth) at 
7 Hz, to remove sky emission; after the first 9-point jiggle, 
the telescope is nodded such that the chop position lies on 
the opposite side of the source. Uranus was used as the pri- 
mary flux calibrator. Skydips were taken regularly to deter- 
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Figure 1. The current status of the high redshift (z >1) 850 fim 
Hubble diagram using the observations presented in Table 2. Ra- 
dio galaxies are shown as filled circles; open circles show the 
submm detected lensed sources with vertical lines connecting the 
observed flux with the unlensed fluxes; open stars show previously 
detected radio quiet z > 4 quasars; filled stars show our new ob- 
servations. Upper limits are plotted at the 2cr limit with arrow 
length equal to the la uncertainty. Also plotted is the expected 
flux from a thermal source like Arp220 as function of redshift for 
two cosmologies using the observed SED from Klaas et al.(1997). 



mine the sky opacity which varied considerably throughout 
the run as reported in Table 1. 

Data reduction was performed with the Starlink surf 
software, and consisted primarily of experimenting with sky- 
removal techniques, using the mean and median of different 
rings of bolometers as a background estimates. The results 
were robust to variations in the precise method, and we take 
the 'mean-of-the-innermost-ring' as our standard sky value. 



3 RESULTS 

At 850 fim, one quasar, BR2237-0607 at z = 4.55, is 
detected with high significance (4.5<r), while the quasar 
BR0019— 1522 has a detection at 2er. Two of the remain- 
ing objects have good lcr w 1.5 mjy upper limits. None of 
the objects are detected at 450 fim; at this wavelength, the 
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Table 2. Physical properties of representative 850 /an sources, assuming Td = 50 K and j3 = 1.5 



Object name Redshift F 850 M d ^FiFtf SFiJ Af s f 

(mjy) (h- 2 10 8 M Q ) (h- 2 10 12 L Q ) (a> X h" 2 M yr" 1 ) 



BR0019-1522 


Quasar 


4.52 


s.sii.e 1 


1.5±0.7 


5.0±2.4 


500 


-27.6 




( 111 V 

Vc^ LI OiJ Cli 


4.52 


< 3.0 1 


/ 1 i 

\ J- .1 


<^ A ^ 


< 450 


—28.2 


p. p. 1 finn-un79Q 


Quels elf 


A ^ 


<- D.Z 


S 1 ^ 
-L.O 




s /ion 


—28 1 


RR9937 nfifl7 


Quels elf 




^ 0—1— 1 1 1 


z.oznu.o 


7 c;-M p. 


7^(1 
t ou 


—28 1 


RRT0Q^9 01 1 c »^- a 
jj jrx.iuyoz — uii J 


v^uasar 


4 43 


14i2 2 


u.oznu.y 


Z ± . Z^O. U 


2120 


97 7 
— z i . / 


BR1033-0327 


Quasar 


4.50 


7±2 2 


3.2±0.9 


10.5±3.0 


1050 


-27.6 


BR1117-1329 


Quasar 


3.96 


13±1 2 


6.4±0.5 


20.8±1.6 


2080 


-28.1 


BR1144-0723 


Quasar 


4.14 


7±2 2 


3.4±1.0 


11.0±3.1 


1100 


-27.5 


BR1202-0725 


Quasar 


4.69 


42±2 2 


19.0±0.9 


62.0±3.0 


6200 


-28.5 


BRI1335-0417 


Quasar 


4.40 


14±1 2 


6.5±0.5 


21.3±1.5 


2130 


-27.3 


8C1435+635 


Radio Galaxy 


4.25 


7.77±0.76 3 


3.7±0.4 


12. Oil. 2 


1200 


-24.4 


4C41.17 


Radio Galaxy 


3.80 


17.4±3.1 4 > m 


7.5±1.3 


24.5±4.4 


2450 


-25.0 


APM08279+5255'' 6 


Quasar 


3.91 


75±4 5 


2.6±0.1 


8.6±0.5 


860 


-31.0 


SMM02399-0136 ! ' C 


Seyfert 2 


2.80 


26±3 6 


6.1±0.7 


19.8±2.3 


1980 


-24.0 


H1413+117 i ' d 


Quasar 


2.55 


66±7 7 ' m 


3.4±0.4 


11.1±1.2 


1110 


-28.5 


F10214+4724'' e 


Seyfert 2 


2.29 


50±5 8 - m 


0.9±0.1 


2.9±0.3 


290 


-23.1 


Arp220 


Starburst 


0.0183 


67000 9 '" 


1.7 


5.5 


550 


-21.3 



Notes: j The uncertainty only reflects the flux uncertainty. | Computed from observed optical continuum magnitudes assuming a 
spectral index of —0.5, except in the cases of the radio galaxies, where i\-band magnitudes were used. 'Lensed objects. Magnification 
factors for the FIR continuum are (a) uncertain; (b) 10( Green & Rowan-Robinson, 1999; see also Downes et al. 1995); (c) 2.5 (Ivison 
et al. , 1998b), (d) 10 (Yun et al., 1997), (e) 30 (Eisenhardt et al. 1996). Derived physical quantities are corrected by these factors. 
m 800/^m fluxes; equivalent 850 (im fluxes would be 15, 56 and 42mjy, respectively. "ISO 170 (im flux (N.B., at z ~ 4, 850 jim 

corresponds to a rest wavelength of ~ 170 /mi). 
References: 1. This work; 2. Buffey et al. (in prep.); 3. Ivison et al. (1998a); 4. Dunlop et al. (1994); 5. Lewis et al. (1998); 6. Ivison et 
al. (1998b); 7. Hughes et al. (1997); 8. Rowan-Robinson et al. (1993); 9. Klaas et al. (1997). 



observations are noisy, because weather was on the whole 
unfavourable, so we disregard data in this band. Our ob- 
served fluxes are compared with those of other cosmological 
submillimetre sources in Figure 1. 

We can derive dust masses and FIR luminosities from 
our observations as described in more detail in McMahon et 
al. (1994). There is considerable uncertainty in these proce- 
dures (e.g. Figure 2), but we can assume physical parameters 
derived for better-studied objects and extrapolate from our 
850 /xm measurements. The mass of a dust cloud optically- 
thin in the FIR is given by 



M d = 



S v (v hs)Dl 



KdBv(Vrcst,Td)(l + z) 



(1) 



where Dl is the luminosity distance and Kd is the absorption 
coefficient: 



= K-S50 X 



/ frest \ ' 

V353 GHz J 



(2) 



We assume a frequency dependence of absorption coefficient 
of P = 1.5, as shown to be a good fit for both low redshift lu- 
minous infrared galaxies (Carico et al. , 1992) and luminous 
high redshift quasars (Benford et al. , in prep.) and ksso = 
0.11m 2 kg -1 , based on the widely-used Hildebrand (1983) 
normalisation at 125 /im(ft=1.875m 2 /kg), which we have ex- 
trapolated to longer rest wavelengths assuming (3=1.5. 

The above assumes that the rest-frame FIR continuum 
is thermal emission from dust which is optically thin to 
FIR photons, which is justified since spectral indices in the 
Rayleigh-Jeans region for typical z > 4 quasars are ^ 3 



(Buffey et al. in prep., Isaak et al. , 1994). It also assumes 
that the dust is isothermal: a single temperature of 50 K is a 
typical value derived from fits to well-covered FIR SEDs. If 
the dust were not isothermal, fits to the currently-available 
data would be biased towards the higher temperatures, caus- 
ing us to underestimat e the true dust m ass. ISO Photome- 
try of nearby ULIRGs ( Klaas et al. 1997 ) suggests that two- 
component fits are more appropriate, low-temperature dust 
at ~50K, presumed to be heated by starburst, and high- 
temperature dust at ~150K, pos sibly heated by an AGN; 
whereas observations of quasars ( |Haas et al. 1998] ) suggest 
much higher (150 K Td 600 K) temperatures for this 
component. Since the existing z > 4 data do not constrain 
two-component fits, we assume that the 850 /im flux is dom- 
inated by the low-temperature component. 

The FIR luminosity is obtained by integrating under a 
single modified blackbody curve. Given, Td=50K and /3=1.5, 
Lfir is directly proportional to 



3.3 x 



M d 



10 s M 



(3) 



Assuming that the FIR flux is due to dust heated by a star- 
burst, the star formation rate is given by: 



SFR = a x 10 



10 ^FIR „ _l 

— — Mq yr 



(4) 



where the value of a depends upon the stellar initial mass 
function (IMF) adopted and the timescale for the burst. 
For example, Scoville & Young (1983) deduce a = 0.77 by 
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Figure 2. Variation of derived dust mass and FIR luminosity 
with dust temperature, relative to values with T = 50 K and (3 = 
1.5, for various values of f3 and absorption coefficient, calculated 
from a 850 ^m flux from a dust cloud at z = 4.5. 



considering the total energy radiated by main sequence O, B 
and A stars [M > 1.6 Mq), over ~ 10 9 yr; whereas Thronson 
& Telesco (1986) obtain a = 2.1, assuming a Salpeter IMF. 
Allowing for the low mass (1.6 M Q > M > 0.1 M ) stars, a 
rises by a further factor ~3. 

Table 2 shows results derived for a sample of cosmo- 
logical submillimetre sources, in a self-consistent manner, to 
enable comparison of observed and derived properties. The 
uncertainties inherent in our assumptions are illustrated in 
Figure 2. Note that the uncertainties in FIR luminosities 
reflect solely the uncertainty in the flux, and do not include 
the uncertainty due to the assumed single temperature. In 
Figure 3 are plotted the inferred FIR luminosities against 
optical luminosity for the z>4 quasars. Since the quasar 
BR1202— 0725 has been resolved into two millimetre and op 



-28.0 

M„ 



1 1 100 

-29.0 
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ticalsources ( Pmont et al. 1996bJ ; |Hu, McMahon & Egami| to Mg = _ 2L0 , the source density reaches ~400 deg 



Figure 3. FIR luminosity and inferred star formation rate 
against absolute B-band magnitude for the z>4 quasars tabu- 
lated in Table 2. Filled squares and upper limits: this work. Open 
squares: previously studied quasars. The dot-dashed horizontal 
line shows the FIR luminosity of Arp220. The dashed line shows 
the locus of Lfir = £tjv- 



We can investigate the AGN contribution to the deep 
SCUBA source counts using a simple model where all AGN 
have a fiducial FIR luminosity of ~ 5 x 10 12 L© with an 
Arp220-like spectrum in the submm range. We assign this 
FIR luminosity to each quasar and adopt a standard Boyle, 
Shanks & Peterson (1988) quasar luminosity function, ex- 
trapolating to high redshifts by assuming no evolution be- 
tween 2 = 2 and z — 3, and a density fall-off by a factor 
of 2 per unit redshift thereafter. Integrating the luminos- 
ity function down to Mb = —23.0, we find that there are 
~140 quasars deg -2 . This absolute magnitude limit is pri- 
marily based on historical convention from ground based im- 
age quality, and if we relax this limit and extrapolate down 

Fi- 



1996) i is plotted as two discrete points, since its extreme 
luminosity might be explained if it represents the merger of 
two gas-rich galaxies or of an AGN and a gas-rich galaxy. 



4 DISCUSSION 

From the results presented in Figure 3, there is no obvious 
correlation between Lfir and Ltjv- Indeed, over a range of 
a factor of 1.5 in UV luminosity, the scatter in Lfir is ~ 10. 
The absence of a strong correlation between Lfir and Luv 
is supported by the results tabulated in Table 2, with the 
obvious caveat that this compilation, is rather heterogeneous 
and consists primarily of detected objects. This conclusion 
has the implication that the high redshift AGN population 
may h ave a FIR luminosity function that is independent of 



nally, if one allows for the observed ratio of Seyferts II to 
Seyferts I at low redshift (2.3T0.7; Huchra & Burg, 1992), 
even neglecting any contribution from totally optically ob- 
scured AGN, it is feasable that that a major part (15-100 
per cent) of the SCUBA source population could possess 
AGN. 

Further evidence for such high space densities of AGN 
comes from the deepest X-ray surveys which reach source 
densities of ~1000 deg" 2 (Hasinger et al. 1998). Whilst the 
nature of the faintest sources is still controversial, an AGN 
origin for the energy source is widely considered as the most 
plausible explanation {e.g. Hasinger et al. 1998). Further- 
more, estimates of the local space density of massive dark ob- 
jects within galaxies (Magorrian et al. 1998) exceed the den- 



sity of the optically-selected quasars, of which these objects 



UV lun Linosity. 



In earli er mm/submm work (tsaak et al. 1994; Omont 
et al. 1996a), only a small fraction (~ 20 percent) of the ob- 



are presumed to be the dull remnants (Haehnelt, Natara- 
jan fc Rees 199£ ; Frentham et al. 1998 ) : there may exist a 



served sources were detectable, and the observational limits 
were at quite high FIR luminosities. These new observations 
indicate that the typical luminous quasar has a FIR lumi- 
nosity less than or equal to that of Arp220. 



significant population of AGN which are either obscured by 
dust or accrete in low-efficiency modes. 

Lower limits on the contribution of AGN to the submm 
source counts comes from studies of local ultraluminous in- 
frared galaxies. These limit the fraction of ULIRGs in which 
AGN activity is the primary energy source to 20-30 percent, 
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although these observations do not exclude the possibility 
that th e remaining starbu rst-dominated cases contain AGN 
as well( |Genzel et al. 199S| ). 

An alternative estimate of the contribution of broad line 
AGN to the submm source counts can be derived by assum- 
ing that Lfir is linearly correlated with Ljjv as shown in 
Figure 3. Above 3mJy at 850 /im, the source density ranges 
between 4 and 9 deg -2 for a thermal spectrum with tem- 
perature of 70K and 30K respectively. This is essentially be- 
cause only AGN with Mb<— 25 would have detectable flux 
at 850 /jm. It is thus important that observations of quasars 
with Mb<— 25 are carried out at 850 /im to distinguish be- 
tween the range of models. 

If AGN are present in a significant fraction of the 
SCUBA sources, this has important implications for studies 
which use the SCUBA source counts to infer star formation 
rates from the from a 850 (jm flux. The tacit assumption in 
Equation 3 is that this locally-derived relation between Lfir. 
and the star formation rate, based on galaxies such as M51 
in which Lfir ~ 10 10 L©, can be applied within the extreme 
environments of high redshift AGN and massive starbursts, 
in which FIR luminosities are 100-1000 times larger. For 
example, the intense UV flux from any starburst or AGN 
could suppress the formation of low mass stars which dom- 
inate the total mass. Thus, even if the dust is not heated 
by a more energy-efficient AGN, the AGN's presence may 
effect the underlying IMF which is used to deduce the total 
star formation rate. 



5 CONCLUSIONS 

We have observed a small but statistically-complete sample 
of luminous high redshift (z > 4) quasars with SCUBA at 
850 and 450 /J,m, and have found that the typical FIR lumi- 
nosity of these quasars is less than or comparable with that 
of the archetypal local ultraluminous IRAS galaxy Arp220. 
Using a simple model for the contribution of AGN to the 
submm source population, we find that i; 15 per cent of the 
sources in the deep SCUBA surveys at the Ss5o~3-4mJy 
level, may contain active nuclei. Hence, both the contribu- 
tion of AGN to the FIR luminosity, and its effect upon the 
underlying IMF (should starburst be the dominant power 
source), need to be considered. 

Further submm studies of lower redshift and lower lu- 
minosity quasars are required in order to determine directly 
the far infrared luminosity function for quasars. In addition 
sensitive X-ray observations with AXAF and XMM may de- 
termine the fraction of the SCUBA sources which harbour 
AGN. 
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